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Abstract 

To leading order in a^, we evaluate the leading and non-leading l/M^ corrections to the B*B7r 
and B*B'y couplings using QCD spectral moment sum rules in the full theory. We find that, for 
large Mb and contrary to the heavy-to-light B — > 7r{p)W form factors, which are dominated by 
the soft light quark vacuum condensate, these couplings are governed by the hard perturbative 
graph, like other heavy-to-heavy transitions. We also find that for the B* B7, the l/M^ cor- 
rection is mainly due to the perturbative and light quark condensate contributions originating 
from the graphs involving the heavy quark part of the electromagnetic current, which are essen- 
tial for explaining the large charge dependence in the observed D*~ — > D~'~f and D*^ —>■ D^j 
decays. Our best numerical predictions without any free parameters for the i?*-meson are: 
dB'-BO-K- — 14±4, F^.-^^^-^ ~ (0.10±0.03) keV and the large charge dependence of the ratio: 
Fst-^B-^/F^.o^^o^ ~ 2.5 . For the Z}*-meson, we find: F£,.-_,£,o^- ^ 1.54F£)*o^£)0^o ^ (8±5) 
keV, F^*-^^-^ ~ (0.09lg:^?) keV and Td*o^do^ ~ (3.7 ± 1.2) keV, where the branching ratios 
agree within the errors with the present data, while the total widths Tiy*o_,aii ~ (11 ± 4) keV 
and TD*-^aii — (12 ± 7) keV are much smaller than the present experimental upper limits. 
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1 Introduction and notations 



The 3*311 and D*Dt[ couplings have been studied by several authors using QCD spectral sum 
rules combined with the soft pion techniques [|I](see also 0), light cone sum rules or heavy 
quark expansion plus soft pion techniques 0, while the 3*3p coupling has been studied recently 
using QCD double exponential sum rules for the three-point function However, though, 
apparently convenient, as one works with the two-point function, the sum rule approach of [Q 
is quite peculiar due to the presence of the unphysical so-called parasitic term, which can only 
be eliminated in a more involved combination of sum rules. Moreover, the connection of the 
light-cone sum rule used in 0, with the standard QCD spectral sum rules (QSSR) involving the 
vacuum condensates, is not crystal clear due to the poor understanding in terms of the vacuum 
condensates of the real structure of the light meson wave functions entering into this approach 
(for a criticism on the unreliability for the construction of the hadronic wave functions on the 
light cone, see e.g. 0). The theoretical evaluations of these couplings are interesting as they 
can be used for determining the normalization of the 3 — * 7r(p) lu form factors near zero pionic 
recoil. Indeed, the QSSR analysis of the g^-dependence of the 3 nW form factor /+ indicates 
that it behaves like a polynomial in g^, which can be fitted to a good approximation with the 
pole form ||^ ^, though one cannot use (like currently done in the literature) such a pole 
parametrization for studying its Mb-behaviour at = Experimental measurements of 
these couplings are expected to be improved and available in the forthcoming high-statistics 3- 
and T-charm-factoTj machines from the processes 3 tttu, D* Dn and 3*{D*) —>■ 3{D)'~f. 
In this paper, we shall use the QSSR double moments sum rule approach in order to study 
the large Mft-behaviour of the previous couplings and to estimate their values. Contrary to the 
popular double exponential sum rule, this approach is quite advantageous in the analysis of the 
three-point function, as it prevents the blow-up of the QCD series when the heavy quark mass 
is large but here the number of derivatives remains finite f^. The couplings are defined as: 

{3*{p)3{p')n{q)) = gs^B.q.e^, {3*{p)3{p'h{q)) = -egs^j^^p^py^^P e,e', , (1) 

where q = p' — p and —Q^ = < 0, while are the polarization of the vector particles. We 
shall be concerned with the vertex function: 

V^''\p,p',q) = -Jd'x d'y e^^'''y-^-\Tfj}^\x)JU^)JB{y)) , (2) 

where the currents are: 

Jl = e<?^7^9 + (^qQi^Q Jl = ("^« + md)u-f^d 

u,d c,b 

J^. = uYb JB = iM, + md)d^'b, (3) 

and u, d, c, b are the quark fields and e^, eg their electric charge in units of e. The vertex 
obeys the double dispersion relation 0: 

vip,p',q)=-^, rT^im^(^'^') +- (4) 

AtT^ JMl S — p^ J A/2 S' — p'^ 

^Here, the dispersion relation is done with respect to the two heavy meson momenta, which is not the 
case of the B tt{p)ID. This different configuration is important for the Mb-behaviour of the different QCD 
contributions. 
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Exploiting the fact that Mi, is much larger than the QCD scale A, where the LHS can be 
evaluated using the Operator Product Expansion a la SVZ |T0[, one can work, in the chiral 
limit, with the double moment sum rule: 

•^'"■"''-i/.^iSr/.;jS^-n.^') (5) 

where n, n' are finite numbers of derivatives evaluated at = p'^ = 0. 



2 Moment sum rule for the B*Bn and D*Dn couplings 

The perturbative QCD expression of the spectral function reads: 

1 T A / X , . ^2 M?(s + s' + Q^) - 2ss' 



4vr2 '^^'"^ ''4vr2^ {(s + ,' + g2)2_4,,.}3/2 

where = — > is the pion momentum squared and where the integration limit condition 
is: 

(s~M^)is'-M^)>Q'Ml (7) 

It is easy to check that, contrary to the case of -B — ttW form factor, where the light quark 
condensate is dominant, the light condensate contribution vanishes here after taking the and 
p'^ derivatives, which is a consequence of the fact that the dispersion relation has been done with 
respect to the heavy quarks momenta like in the case of a heavy-to-heavy transition. The other 
remaining effects which are suppressed by 1/ compared to the leading perturbative diagram 
will be neglected to the approximation we are working ^ The phenomenological side of the 
sum rule is parametrized using the usual duality ansatz: lowest resonance +QCD continuum 
from the thresholds Sc and s'^. By transferring this QCD continuum effect into the QCD part 
of the sum rule, one obtains 0: 

,An,n')_ V2MB*fB* v/2M|/b y2m^^ 1 ds K ds' 

Jd* Jd " 

where the coupling constants are normalized as: 

(0| JlK) = V2UM'^ {0\J^\p) = V2-^e^, 

(0|Jb|5) = v^/bM| (0|J^.|5*) = v^/B*MB.e^ (9) 

where = 93.3 MeV and 7p = 2.56. In the case where oo (static limit), it is convenient 

to work with the non-relativistic variables E and E' defined as: 

s = {E + Mb)^ and s' = {E' + M^f , (10) 



^In this region of > 0, the question of non-Landau and complex singularities and of anomalous thresholds 



do not arise |ll . 

^One can however notice that the four-quark condensate contribution behaves like l/Q"* which reflects the 
fact that the present approach cannot be used at = as expected. 

^Here and in the following, we shall neglect the contribution of the 7r'(1.3) similarly to previous analysis of 
the ujpTT- and TrA'^iV-couplings using vertex sum rules pM- 
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and to introduce the new variables: 

x = E-E' and y = ^{E + E') . (11) 

Due to the almost good symmetry between the B and the B*, we shall use: 

Mb* ~ Mb, E, ~ E^, n = n' = . (12) 
By keeping the non-leading l/M^-terms in the expansion, we obtain to leading order in a^: 

Mi"'""' ~ (mu + md) \u ^Q / . 



^4(n3+l)^2^ 7q Vi^^W (a;2 + Q2)3/2 

^ ' - ^{Q' + + 2l/^) - 2(2n3 + l)y + ^5(y - i^^^ + Q^) 



Mr 



■ (13) 



For consistency, we shall use in our analysis the lowest order expression in as of the decay 
constants from the moment sum rules [|l^]: 

- 2Tr^MB\Mj Y 2^''' + ^^M, 2 E^ ] 
fl - + (14) 



27t'^ Mb \ Mb J [ 2' 3' Mb 2 _ 

consistent with the normalization of the currents in Eq.(3) and with the definitions in Eq.(9). 
One should notice that the overall {Mb /Mb) factor also brings a l/M^-correction which tends 
to reduce the apparently huge correction in the curly brackets and leads after the moment 
sum rules analysis of Jb to the well-known 1 GeV/Mb-correction to this quantity. One can 
also notice that the l/Mf,-correction to J'b* is slightly smaller than the one of Jb as generally 
expected. Using the previous formulae in Eq. (14), the emerging effective values of Ec fixed 
from the numerical analysis of Jb and Jb* including the corrections are [l^, p!3| : 



~ (1.6 ± 0.1) GeV ~ (1.3±0.1) GeV ~ (1.1 ± 0.2) GeV , (15) 

which, using Eq. (10), can be parametrized as: 

^-^~(^-2|)- <«' 

As in P], we minimize the n-dependence of the results by requiring that the leading term is 
n-independent. This leads to the constraint: 

4^3 + 1 = 2^2 + 1 , (17) 

where n2 — 4 — 5 is the value where Jb from the two-point function has been optimized [0 . 
By evaluating numerically the different integrals, we obtain to a good approximation: 

mi 3 i?c iuu) 1 , , 



where: 



g'^B*B. - ^7 -^^Mb <Io = i^ dx ^ . (19) 



The analytic expression of the integral Xq is: 



Up = Q/Ec) = |p 



3 ^2 



1 + |P 



- p 



(20) 



which exhibits a broad maximum in the range 1-3 GeV^. At — 2 GeV^, where the ab- 
solute maximum is obtained, its numerical value for different can be parametrized by the 
interpolating formula: 

Jo ^ (0.119 ±0.001)£;c ■ (21) 
Therefore, using Eq. (16), we finally obtain: Q 

where we have introduced the static coupling g°°: 

(0.119^;?°) , (23) 




which controls the interaction of the pion with infinitely heavy fields in the effective Lagrangian 
approach: 

Ant = \g^TrH^^^,{n^d,n - 7rd,n^)H, (24) 

where H and tt are the heavy and pion fields. The M;,-behaviour obtained here, which is 
dictated by the one of is in agreement with current expectations Il|-0. The agreement with 
the one in [jl[ (the sum rule used in is very similar to the light-cone sum rule in the treatment 
of the pion) and the light-cone sum rule [Q] can be mainly due to the fact that, in the present 
process, the hard perturbative diagram gives the leading contribution in 1/Mh, where the 
present version of the light-cone sum rule approach, which is dominated, (by construction), by 
the hard perturbative diagram, is appropriate. In the case where the soft process is dominant, 
like e.g. in the analysis of the B 7i{p) semi-leptonic processes [§, one should need a modified 
version of the light-cone sum rule approach in order to take properly into account the dominant 
non-perturbative (uu) condensate contribution. We expect that, in the moment sum rules , the 
as correction is much smaller than the one in the non-relativistic exponential one, as here is 
evaluated at a larger scale of about MfeZ-^n, while in the exponential sum rule the scale is much 
lower at about 1 GeV. Therefore, we expect that the expression in Eq. (22) with the value of 
Ec in Eq. (15) gives a good approximation of the physical result. However, we consider, as 
an intrinsic error of the approach, the known 30% effect due to as in the decay constants from 



the moment sum rules |T^. We shall use (m„ + m^^) (IGeV) =(12.5±2.5) MeV [0 rescaled at 



= 2 GeV^, and the corresponding value of the quark condensate. Then, we deduce: 

~ (0.15 ±0.03). (25) 

where the error takes into account the effect of 30% by the radiative correction to the value of 
J'b [0. Our prediction in Eq. (25) is in agreement with the range of values obtained in |l|. 



though the authors in use a too low value of f^"'^^'^ ~ 1.45/7r, while we use here the two-loop 

^We have checked using the complete perturbative expressions of the three- and two-point functions that the 
higher order 1 /M corrections are small and do not spoil the validity of the following approximate formula even 
at the charm mass. 
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value ~ 2/,^ from and from the recent lattice results ||15[. For the physical 5- meson, 



and including the 1/Mi, correction which is of the order of +28 % at the 6-mass, we obtain: 

gB*Bn ^ 14.5 ± 3.3 , (26) 

where again the error takes into account the effect of radiative corrections to Jb- We have used 
the two- loop non-relativistic pole masses ||16||: 



Mb = (4.7 ± 0.03) GeV M, = (1.45 ± 0.05) GeV, (27) 

consistent with the present use of non-relativistic sum rules. One can notice that, like in the 
case of /b, the l/M;,-correction is large (28% at the 6-quark mass and 76% at the c-quark mass). 
This feature can make the extrapolation of the result to the D-meson quite risky. However, 
as already mentioned earlier, our explicit evaluation of the complete perturbative expression 
of the three- and two-point correlators indicate that higher order corrections in 1/M remain 
small. Therefore, we can deduce, with a quite good confidence, the estimate from the moments: 

gD*Dn ^ 7.1 ± 1.6 . (28) 

We cross check the validity of the previous result by invoking semi- local duality sum rules for 
the two- or three-point functions which correspond respectively to the particular cases where 
(n3,n2) = (—1/2, —1) or (^3,722) = (—1, —2), and which have been discussed extensively in the 
case of the QCD two-point functions for light [|l8l and heavy quarks |jl9|, |T^. In these particular 



cases, the 1/M corrections to the leading term of the three-point function are much smaller 
and result by a correction of about +(9-ll)% and -|- (23-26)% for the coupling respectively at 
the 6 and c quark masses, giving []: 

gB'B. ^12.7 ±2.9 , gD'Dn^5.0±l.l . (29) 

By considering the previous results, we conclude that the most conservative estimate of the 
couplings from the sum rule is: 

~14±4, go^Dn - 6-3 ± 1.9 . (30) 

Our final results are in good agreement with the ones in |jl| but are much smaller than the 
ones from and to the indirect determination of which is correlated to a higher input 



value of the B irlu form factor. Taking into account that the l/M^-correction to is 
approximately equal in strength but opposite in sign with the one for gs's-K, one obtains with 
a good accuracy: 

„ _ gB*BTTfBWMB , s 

Rbd = 7 — 7fr= - 1 > (31) 

as expected from an alternative analysis [Q. Our numerical values of the physical couplings 
are in better agreement with the results in [jl| including the radiative corrections than with the 
ones in [^, which is higher than ours by a factor 2 . f\ The coupling in Eq. (30) leads to the 
prediction: 

Td'-^do^- = JnS-l^-l' - 1.54r^*o_^.vo ^ (8 ± 5) keV , (32) 



^One should also notice that the use of the Laplace sum rules leads to a small value of the optimization scale 
r, which is pratically similar to the semi-local duality sum rules used here. 

*The agreement with Q for the physical B-meson is due to the approximately same value of fs used here 
(I.S/tt) and in (1.36/^), though we, orginally, do not start with the same value of In our analysis, 

the value of /^''^* ~ I-S/tt, at the physical i3- meson mass has been deduced from /^*°* ~ 2/,^, after taking into 
account the l/Aff,-correction. In Ref. this l/Aff,-effect seems to be much smaller than currently expected 
as, there, f§-' ~ 1.45/. « fj^^'. 
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where we have assumed isospin invariance for the couphngs. Using the observed branching 
ratios [0], one can also predict the radiative decays: 

r^.o^^o^ - (3.0 ± 1.2) keV TD,-_^D-y ^ -l^-H keV . (33) 

and the total widths: 

r^*-^,^^ ^ (12 ± 7) keV r^.o^,;; ~ (8 ± 3) keV . (34) 

The predictions for the total widths are much smaller than the present experimental upper 
limits. An improved measurement of the Z)*-total widths in the next tau-charm factory machine 
should provide a decisive test for the validity of these extrapolated predictions. 



3 Moment sum rule for the B*B^ and D*D^ couplings 

The QCD expression can be decomposed into a light (g = m, (i, s) and heavy {Q = c, h) 
quark parts. For the corresponding vertex function, the one available in ^ agrees with our 
recomputation apart for the relative sign between the perturbative and quark condensate con- 
tributions in the heavy quark component of the electromagnetic cuurent. After a systematic 
l/Mft-expansion of the full QCD expression, one can inspect that the dominant contribution 
comes from the perturbative graph related to the light quarks coupled to the electromagnetic 
current. The heavy quark contribution is l/M^-suppressed compared to the light quark one. 
However, the perturbative and light quark condensate contributions are of the same order in 
l/Mft in this heavy quark component. By keeping the 1/Mb-correction, the QCD part of the 
sum rule reads: 



QCD - ^^M^n^+l) ^^^\2^ 



dx 



y dy 



i^/^W (x2 + Q 



2^1 3/2 



16 asjuu)'^ eg 
9 m5"^+') 



2(2^3 + l)y + 



1 + 
1 



{uu) + Oil/Mb) 



(35) 



where eg(Q) is the charge of the light (heavy) quark in units of e. The phenomenological side 
of the sum rule can be parametrized as: 



phen 



gB-B-y- 



V2MB*fB' V2M|/b 



M; 



2(n3+l) 



B 



Using an approach similar to the one done for B*B7i, we deduce the sum rule: 

gB*B,{Q')^9i^B,iQ')+9B^B,iQ') 



(36) 



(37) 



where: 



9'B*B,iQ') 
9B*B'yiQ^) 



eg 
Mb 



3 



, 2Mb, 



16 as{uu)'' 
-n- 



9 {E^rQ' 



B\^ 



(38) 
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where Xq is the integral defined in Eq. (19). For oo, the couphng is given by the hght 

quarks contribution and remains constant. The l/M^-correction in the hght quark contribution 
is much smaller than the one for gB*Bn since there is an almost cancellation of the l/Mf,- 
correction with the one from the £'c-dependence of Xq as can be deduced from Eq. (16), while 
the one due to the heavy quark is important at the c-quark mass. The light quark coupling 
exhibits a typical monopole behaviour for > , while the heavy quark coupling is Q^- 
independent. Therefore, we use a light vector meson dominance for the estimate of the light 
quark coupling, which can be related to the 5*_BV^-coupling (V = p, u) as: 

M') - (^) . (39) 

A sum rule analysis of the S* 5 V^- coupling similar to the one for B*Btt, shows a very good 
stability for 0.4 < < 2.2 GeV^, where the optimal value obtained for ^ 1 GeV^ reads: 



'V2MI 
27v 



^^i^v- ^ (0.84 ± 0.10)/ (40) 



We have used as{uu)'^ ~ (5.8 ± 0.9)10 ^ GeV^ ||2T|, [T^, which shows a negligible contribution 
of the four-quark condensate in the range of Q^-stability. Therefore, we deduce: 



gB'B^iQ^ = 0) 
gD*D^iQ^ = 0) 



, n .rN (0.90 ±0.16) GeV 
e,(1.14 ± 0.15) + cq'- 

(0.90 ±0.16) GeV 
e,(l.ll±0.24)±eQ^ ' 



GeV"^ , 

GeV"^ . (41) 



For the 5-meson, the heavy quark contribution is relatively small. One obtains: 

Tb*--.b-, = gl^Bj^W-yf - (-10 ± .03) keV . (42) 
and with a better accuracy for the ratio: 



B*o-,BO'r 



2.5 , (43) 



which deviates strongly from the naive static limit (M^ oo) expectation (e„/erf)^ = 4. For 
the D-meson, the heavy quark contribution is relatively important. One should notice that the 
T-contribution has been completely ignored in the phenomenological analysis of , which can 
explain their opposite prediction of this ratio with respect to us and to the data. We deduce 
within the previous approximations: 

Td*o^dOj ^ (8.0 ± 2.7) keV r^._^^_^ ~ (.01 ± 0.08) keV . (44) 

If one instead works with the semi-local duality like-sum rule using the complete expressions 
of the perturbative contributions, one finds that the heavy quark contribution to the coupling 
is reduced by 30% and leads to: 

Td*o^dOj ^ (6.7 ± 2.4) keV 1^*-^^-^ ~ (.04 ± 0.08) keV . (45) 

Therefore, the most conservative sum rule estimate is: 

Td*o^dOj ^ (7.3 ± 2.7) keV Td*-^d-^ - (.03 ± 0.08) keV . (46) 
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which, despite the large error, shows that the heavy quark contribution acts in the right direction 



for explaining the large charge dependence of the observed decay rates [1/]. One can combine 



these results with the ones in Eq. (32), for an attempt to deduce the ratios of rates: 

r£)*o^2)07ro/rz).o^£)o^ r£)*-^£)o^_/r£i.-^£)-^, (47) 

but, due to the large errors, the comparison of the predictions with the data is not very con- 
clusive. Alternatively, we can combine the predictions in Eq. (46) with the observed branching 
ratios given in Then, we predict: 



r^.o^flo^o ~ (14 ± 5) keV Tn*-^DO^- < 18 keV. (48) 

and: 

r^.o^,,, ~ (20 ± 7) keV rn*-^aii < 27 keV . (49) 

These results are respectively in fair agreement within the errors with the direct calculation in 
Eq. (32) and with the prediction in Eq. (34), though the ones for the D*~ have a large error due 
to the inaccuracy of the measured and predicted D*~ branching ratio. The agreements 

between the different results given in this paper is an indication for the self-consistency of the 
whole approach. 



Conclusion 



We have systematically studied the couplings P*P7t and P*P'y (P = B,D) using a l/Mh- 
expansion in full QCD with the help of moments sum rules. It is important to notice that, 
like other heavy-to-heavy transitions, the couplings are dominated by the hard perturbative 
diagram. This is not the case of the B 7r(p) lu and B K*'y heavy-to-light transitions 
which are governed by the soft light quark vacuum condensate @]. Technically, this difference 
is mainly due to the uses of different dispersion variables for the heavy-to-heavy and heavy- 
to-light transition processes. We find that, for the P*P'k couplings, the l/Mf,-corrections due 
mainly to the perturbative graph are large but they tend to cancel for the quantity fpgp*pn and 
implies to a good approximation the relation in Eq. (31). For the P*P7-coupling, the 1/Mb- 
correction is due mainly to the perturbative and light quark condensate contributions from the 
heavy quark component of the electromagnetic current, which goes in the good direction for 
explaining the large charge dependence of the ratio of the — > D^'j over the D*~ —>■ Z)~7 
observed widths. For the P*-mesons, our predictions are given in Eqs. (30), (42) and (43), 
where for experimental interests in the next P-factory machine: 

Tb''~-.b-^ ^ (.10±.03) keV , 

TB*-^B-^/'^B*o^BO-y - 2.5 , (50) 

where the latter deviates strongly from the naive static limit {Mf, oo) expectation (e„/erf)^ = 
4. By combining the previous different results of the P*-meson, our averaged predictions for 
the different exclusive widths are: 

r^.-^^o^- ~ 1.54rB,o^Bo^o ~ (8 ± 5) keV. (51) 

and: 

TD*o^DO-y ^ (3.7 ± 1.2) keV TD^-^D-y - i-OQtfor) keV . 

TD'0_,,a ~ (ll±4)keV rz?*-.^,^ ~ (12 ± 7) keV . (52) 
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The branching ratios agree within the errors with the present data though the total widths are 



well below the experimental upper limits F/j.-^a// < 131 keV and Tj:,*o^aii 

< 2 MeV [Ol. We 

urge experimentalists to improve the measurements of these total widths in the near future, 
as these measurements are necessary for clarifying the present disagreements between different 
theoretical predictions. 
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